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Abstract 

Within an early Universe scenario, nonextensive thermostatistics (NET) is 
investigated on the basis of data concerning primordial Helium abundance. 
We obtain first order corrections to the energy densities and weak interactions 
rates, and use them to compute the deviation in the primordial Helium abun- 
dance. After compare with observational results, a severe bound is stablished. 
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For a variety of physical reasons much work is nowadays devoted to nonlinear formalisms. 
Among them we can single out nonextensive thermostatistics (NET) [|TJ as an active area of 
research. NET is based upon the following two postulates @.||]: 

Postulate 1.- The entropy of a system that can be found with probability pi in any of 
W different microstates i is given by 

w 

s g = (q - 1)- 1 Up* ~v% (1) 

i=l 

with q a real parameter. We have a different statistics for every possible q— value. Of course, 

5> = 1 > ( 2 ) 

i 

and it is easy to see that, for q = 1, one regains the Boltzmann-Gibbs form 0. The resulting 
physics is extensive just for q — 1. Otherwise we are led into the realm of nonextensivity 

Postulate 2.- An experimental measurement of an observable A, whose expectation value 

in microstate i is ai, yields the q- expectation value (generalized expectation value (GEV)) 

w 

< A > q =J2p 9 i a i , (3) 
i=i 

for the observable A. 

Unappealing as the above postulates may perhaps be considered, it should be strongly 
stressed that these two statements have the rank of axioms. As such, their validity is to be 
decided exclusively not by vague discomfort feelings but by the comparison with experiment 
of the conclusions to which they lead. One such a test was recently reported in refs. 
where bounds to \q — 1| were stablished using the cosmic blackbody radiation. Also, present 
day determination of Stephan-Boltzmann constant a puts a similar constraint, for which the 
order of magnitud is \q — 1| < 10~ 4 . However, it was later noted J7|, that in both cases, the 
application of thermodynamics to these contexts is strictly local, and thus, a non-violation of 
non-extensivity in a large scale could not be sustained on this basis. In this communication, 
we intend to find bounds on non-extensivity not affected by such criticism, which we deem 
fair. 

2 



The phenomenal success of thermodynamics and statistical physics crucially depends 
upon certain mathematical relationships involving energy and entropy, and much work has 
been devoted to i) show that many of these relationships are valid for arbitrary q, and ii) 
to find appropriate generalizations for the rest. In this vein we just mention that, by suit- 
ably maximizing ([[]), Curado and Tsallis || found that the whole mathematical (Legendre- 
transform based) structure of thermodynamics becomes invariant under a change of the 
q— value (from unity to any other real number), while the connection of NET both with 
quantum mechanics and with Information Theory was established in [|J, where it was shown 
that all of the conventional Jaynes-Boltzmann-Gibbs || results generalize to the Tsallis' en- 
vironment. For more details see []]]. Of course, to verify that NET is useful, it is necessary to 
show that it appropriately describes certain physical systems with q- values that are different 
from unity. Much work in this respect has been been performed recently. We may cite 
applications to astrophysical problems P,[10|, to Levi flights [|ll|], to turbulence phenomena 



|12|| , to simulated annealing etc. The interested reader is referred to [|I|] for additional 



references. Now, NET establishes a different (from the orthodox) fashion of doing 'statis- 
tics', i.e., a non-conventional way of counting, that has proved to be useful in a variety of 
contexts. The difference is governed by the value of the Tsallis parameter q. It is clearly 
recommended to use data concerning diverse natural phenomena to estimate the g-value 
with reference to different scenarios. In the present communication we shall attempt to use 
early Universe Helium abundance data so as to find one of such estimates. 
For the canonical ensemble, (P gives [|J 
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1 - (1 - q)(3H\ x -« , (4) 
as the appropriate density operator and 



Z q = Tr 



i - (1 - q)(3H] ^ , (5) 



as the associated generalized partition function ||. Here, as usual, (3 = 1/kT and H is the 
hamiltonian of the system. 



We shall focus our attention upon the f3(q — 1) — + limit, in which a first order expansion 
allows for analytical computations. The expression of the generalized mean value of an 
operator was computed in this limit by Tsallis, Sa Barreto and Loh 0. When applied to 
particle number operators, the concomitant result reads 



< n > g =< h > BG Z%q 



1 + (1 - q)x 



< n 2 >bg , x f 2 <n 3 > BG 

77 < n >BG 



(6) 



< n > BG 2 \ ' <n > BG / 

where x stands for e/kT (e is the energy of a single particle) and the symbol BG means to 
be computed within Boltzmann-Gibbs' statistical tenets. 

With the standard values of < n 2 >bg an d < n 3 >bg, for fermions and bosons, the 
corrections to the energy density in the early universe may be computed fL4| . When the 
particles are higly relativistic, T>m, and non-degenerate T ^> /i, we get 



Pf 



ermions 



2ir 2 Jo 
2tt 2 Jo 



POO 

/ dEE 3 < h > bosons ^ (7) 

J 



POO 

JO dEE * <fl> fe~s, q , (8) 



gbj stands for the degeneracy factor of each one of the species involved. Using (|5]), we finally 
obtain 

7T 2 1 

Ptotal = Pbosons + P fermions = 1^9^ + — ^ (40.02^ + 34.70^/) T 4 (<? - 1), (9) 

3U Z7T 

where g = gj, + 7/8(7/. At high enough temperatures, the energy density of the universe 
is essentially dominated by e~,e + ,u and V. Interactions among these particles keep all of 
them at nearly the same temperature. Accordingly, we set g b = 2 and gj = 2 + 2 + 2x3 
and reach thus the final form of the Tsallis correction to the energy density, namely 

Ptotal = Pstandard + 21.63 T 4 (g - 1). (10) 

We turn now our attention to the details of the computation of those corrections due 
to the weak interaction rate. This rate allows one to compute the neutron abundance as 
the universe evolves. We shall denote by X pn (T) the rate for the weak processes to convert 



protons into neutrons and by X np (T) the rate for the associated, reverse ones. Following the 
standard computations |]I5|,|1S|| , it is possible to see that, for high enough temperatures, the 
weak interaction rate is A(T) ~ X np + X pn ~ G 2 F T 5 , with X np being related with X pn by the 
principle of detailed balance JTB[: X np = exp (— Q/T) X pn , Q = m n — m p = 1.29MeV and 
Gp the Fermi constant. We want to compute first order Tsallis corrections ((q — l)-order) 
to A(T). To do this, we need to analyse the individual interaction rates, i.e., each one of 
the terms in the sum 



Xnp ^u+n—^p+e^ 



-v + X 



(11) 



given by fL5f , 



A 



v+n— >p+e~ 



x f 



A / dp u plp e E e (l- < n e >) < h v > 
Jo 



roo 

, p+0 = A dpeP%p u E u (l- < h v >) < h e > 
Jo 



X 



n—tp+e +v 



fPO 

A / dp e p e p u E u (l- < h v >)(1- < h e >) 
Jo 



(12) 



(13) 



(14) 



where A is a constant fixed by the experimental value of X n ^ p+e - +u . In the preceding 
equations, we have to consider, of course, energy conservation as a boundary condition 
which relates E v and E e . Due to the high temperature limit we are interested in, several 
approximations are in order 

i) All temperatures involved in the present game will be taken as equal, T e = T 1 = T v = 
T, which ensures that reverse reactions have the same form as the direct ones. 

ii) We shall neglect Pauli factors (1— < h v >) and (1— < n e >), and also the electron 
mass in (HI and (IT3|). 



With these approximations (|i~2"|) and ( |i~3"D become identical. Using p e = E e = Q + E u in 
(|T2|), the standard result follows. 

Passing now to the nonextensive context, we must consistently use the < h > q distribu- 
tions functions. Performing the previous integrations, we obtain the leading order corrections 
terms in the fashion 



A, +rwp+e - = KTnXte- + (480T 5 + 2 x 84T 4 Q + 18T 3 Q 2 ) (1 - q)A (15) 

where, 

K+n%+e- = ( 4!t2 + 2 x Z\TQ + 2!Q 2 ) AT 3 . (16) 

In order to get some fresh insight into the problem, we shall consider here only the 
first correction, proportional to T 5 . A more detailed analysis of these weak rates and the 
neutron-proton abundance ratio they yield is considered elsewhere [17]. As explained in [16 



the high temperature regime makes X np ~ 2 x \ u+n ^ p+e - and A ~ 2A np . As a consequence, 
the change in the weak reaction rate adopts the form 

^ = 1920 T 5 (1-q). (17) 

We have now all the ingredients of the nucleosynthesis recipe at our disposal. Basically, 
nucleosynthesis is the competition between the weak interaction rate and the expansion 
rate, given by the Hubble constant via the Einstein equations. The A He production may be 
estimated -in the standard Big Bang model- as 

Y - = A (if;),, < 18 > 

where A = exp(— (t nuc — tf)/r) stands for the fraction of neutrons which decayed into pro- 
tons between tf and t nuc , with tf (t nuc ) the time of freeze out of the weak interactions 
(nucleosynthesis), r the neutron mean lifetime, and x = exp(—Q/kT) the neutron to proton 
equilibrium ratio |L8|]. It is straightforward to compute the deviation produced in Y p by a 
variation in Tf, and correspondingly, tf. We get 

" 6T f 



SY p = Y p 



2XJ \Y P 



(19) 



where a radiation era relationship between time and temperature of the form (T oc t~z) 
is assumed and the one puts 5T nuc = 0, because it is fixed by the binding energy of the 
deuteron. Similar studies concerning bounds on gravitational theories were analized in 
fL9^f2T|. Considering now, Y p = Y° hs = 0.23 and SY P = 0.01, which is the observational error 
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22|, and standard values for the times and the mean life of neutron -which in fact, is not 



modified at order (q — 1)-, we must ask for 

ST f 

0.01 > 0.3766|— ±\ (20) 
T f 

to be satisfied in order to get an estimate of primordial 4 He production compatible with 
observational data. 

In order to obtain a value for -jr- in nonextensive statistics, we equate 



fa\ 8itG 

A K (J = V— (21) 

Adding up the corrections due to the changes in the energy density and in the weak 
interaction rate, our first order result up to (q — 1) reads 

= 1276.4( 9 -1), (22) 

which allows for a stringent bound, using (pC|) , on the value of q 

\q-l\ < 2.08 x 10" 5 . (23) 

In obtaining this bound we reach the main goal of the present communication: the 
early universe physics places a bound upon the Tsallis parameter q. Since the measured 
value of Y p comes from a sample which has been thoroughly mixed, at least during the 
life of the galaxy, and our estimate of q is based on a volume of the order of the horizon at 
nucleosynthesis epoch, the present test avoides the locality problem pointed out in ref . @ . Of 
course, for other kinds of processes, different bounds may be obtained. Our result pertains 
to those taking place during the early childhood of our Universe. Thus, non-extensivity 
is severely constrained upon all epochs of cosmic evolution, with separate, independently 
tested, observational evidence. 

The authors acknoledge partial support from CONICET and valuable conversations with 
D.E. Barraco. 
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